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We study the generation of non-equilibrium spin currents in systems with spatially-inhomogeneous 
magnetic potentials. For sufficiently high current densities, the spin polarization can be transported 
over distances significantly exceeding the intrinsic spin-diffusion length. This enables applications 
that are impossible within the conventional spin-diffusion regime. Specifically, we propose dc mea- 
surement schemes for the carrier spin relaxation times, Ti and T2, as well as demonstrate the 
possibility of spin species separation by driving current through a region with an inhomogeneous 
magnetic potential. 



The rapidly developing field of spin-sensitive electron- 
ics promises to reinvigorate conventional semiconductor 
electronics [|]. Utilizing the carriers' spin degrees of 
freedom in addition to their charge opens possibilities 
for conceptually new spintronic devices. Already im- 
plemented are GMR read-heads and non-volatile mag- 
netic RAM (MRAM); there are also many proposals for 
other devices |2j. Much of the promise of spintronic de- 
vices lies in the enhanced flexibility of the spin degrees 
of freedom compared to electrical charge. While creating 
local charge imbalance generally involves large electro- 
static Coulomb energies, non-equilibrium spin polariza- 
tions can be created at a very low cost. As a consequence, 
spin polarizations are inherently susceptible to manipula- 
tion, e.g. by external magnetic fields or by current drive. 
Indeed, it has been observed experimentally that non- 
equilibrium spin polarization can be transported over dis- 
tances exceeding 100 micrometers with only moder- 
ate diffusive spread. This would be impossible in the 
case of the non-equilibrium charge inhomogeneity, since 
the Coulomb repulsion would blow such inhomogeneity 
apart almost instantly. 

The immunity of the spin polarization with respect 
to electrical screening makes it possible to satisfy drift 
transport conditions, for which the local spin polariza- 
tion can be carried downstream by an electrical current 
over distances proportional to the drift velocity. This 
situation can be obtained for sufficiently large current 
densities, when the spin-drift length becomes larger than 
the spin- diffusion length (both defined below). In this 
Letter we explore this new regime, readily achievable in 
spin transport, but almost never realized in electron-hole 
transport. It is studied based on general formulation for 
spin transport using spin drift-diffusion equations. 

The spin-drift regime is generally favorable for spin- 
tronic purposes as it can provide larger spin-injection 
depths compared to the spin-diffusive regime. Here we 
propose two specific applications of the spin-drift regime: 
(1) dc transport-based measurement of longitudinal (Ti) 
and transverse (T2) spin relaxation times; (2) spin species 
separation by inhomogeneous magnetic fields. To deter- 
mine the spin relaxation times, we propose to measure 



the spin polarization decay profile in the presence of con- 
stant current; the decay length can then be related ei- 
ther to Ti or T2, depending on the experimental setup. 
Unlike conventional pump-probe techniques, there is no 
need for ultrafast time-resolved detection, and the spin 
polarization profile can be measured either by optical |^ 
or some other techniques. To spatially separate spin- 
up from spin-down carriers, we propose to drive current 
in a 2D geometry containing a region of inhomogeneous 
magnetic field. This technique can be used as a flexible 
alternative to the standard spin injection techniques. 

To study spin transport we use the drift-diffusion ap- 
proximation for spin currents in the presence of an inho- 
mogeneous spin-dependent potential Va{v) — ~gfiB<y 'H- 
At this point we assume that in the plane of the elec- 
tron gas, the applied magnetic field H can vary in value 
but not in orientation. Then all carriers in the sample 
can be divided into spin-up (cr = 1/2) and spin-down 
(cr = —1/2) species that feel magnetic potentials equal in 
magnitude but opposite in sign, U|(r) = — V|(r). For the 
external potential smaller than the thermal energy in the 
non-degenerate case, Vcr{r) < ksT, or less than the Fermi 
energy in the degenerate case, Va (r) < Ep , the local car- 
rier density will remain essentially unchanged, both with 
and without current drive. Hence, in what follows we 
neglect electrical screening effects since they only mod- 
ify the results in the order {Va /max{kBT, Ep))'^ . In the 
drift-diffusion approximation, the spin-polarized current 
is 



(1) 



where q = ±e is the carrier charge, the number den- 
sity of the spin species cr, D the diffusion constant, fi 
the mobility, E the external electric field, and V denotes 
the gradient operator. We assumed that the transport is 
unipolar (e.g. n-doped semiconductor), and the diffusion 
constant is the same for both spin species 0| . We do not 
include the Lorentz force, which also acts on the mov- 
ing electrons, since it affects both spin species identically 
and, moreover, is compensated through the classical Hall 
effect for a confined sample geometry. In practice, lo- 
cal variation of the magnetic potential can be achieved 
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either through the use of micromagnets that can create 
strongly inhomogeneous external magnetic fields, or by 
magnetically doping semiconductors, which can produce 
strong position dependence in the carrier g-factor 

The dynamically generated non-equilibrium spin po- 
larization relaxes towards its equilibrium value. Under 
constant drive, the relaxing spins are replenished by the 
divergence of the spin current. In the relaxation-time 
approximation, 



v-j. = -(5,/ri, 

ria = nQex.p{-Va/kBT)/2. 



(2) 
(3) 
(4) 



The expression for the equilibrium density, ria-, corre- 
sponds to non-degenerate carriers. For degenerate car- 
riers, the Boltzman distribution should be replaced by 
the Fermi-Dirac distribution. In the non-stationary case, 
there is an extra time derivative, dSa/dt, on the left hand 
side of the continuity equation. Combining Eqs. (|l|) and 
, the dc equation for spin density under uniform drive 
(V • E = 0) is 



charge imbalance, and hence is screened out for distances 
exceeding the Coulomb screening length . The solution 
of the homogeneous equation for is 



Scr{x) — Aexp(Aia;) -t- i3exp(A2a; 

Al(2) 



(7) 
(8) 



Matching the (x > 0) and (x < 0) solutions at x = 
and requiring vanishing So- at infinity, the spin density, 
s = — Si, for the positive drift velocity {qE > 0) is 



S 



^ exp(Aix), X < 0, 
+ I exp(A2x), X > • 

The solution for the negative drift velocity can be ob- 
tained from this solution as — s(— x). 

Two interesting special cases correspond to the spin- 
drift [Le S> Ld) and the spin- diffusion {Le < Le) 
regimes. In the spin-drift case. 



q 



DV'S, = (5) 



The first term on the left hand side represents the drag 
of the spin density by the electric field (/iE = Vd - the 
drift velocity), the second term, as we will see, is usu- 
ally unimportant, and the last term represents diffusion. 
Similar descriptions have been previously considered by 
several authors ||^, however, with the focus on the dif- 
fusive regime (neglecting E(5cr terms). The drift terms 
have only been invoked in the spatially localized regions, 
such as depletion layers and Schottky barriers Q. The 
primary goal of the present work is to explore the spin- 
drift regime that occurs for larger drives, that is when the 
spin-drift length, Le = jiETi, exceeds the spin diffusion 
length, Ld = VDT\. 

We first consider the one dimensional case with spin- 
polarized negative half-space. The static equilibrium 
density is 



no/2 = aSQ{ 



(6) 



where Q is the step function equal to unity for positive 
arguments and zero otherwise; S is the static spin density 
in the polarized region. The step in the density profile 
implies a step in the magnetic potential, Va- The density 
modification in the presence of a uniform current drive 
is described by the ID version of Eq. (|). The strongest 
singularity is caused by the jump in Vo- at x = (sec- 
ond term in Eq. (H)), which implies discontinuity in the 
density, S^, 



64+0) - 64-0) _ n 



64+0) + 64-0) 2qD 



(l^,(-0)-l-,(+0)). 



Within our assumptions about the magnitude of the 
magnetic potential, Va < max{kBT, Ep), this correc- 
tion is small and hence can be neglected. Also, it causes 



5 exp(— x/L^), X > 0, 
0, X < 



(10) 



And in the spin-diffusion case 
SLe 



2L 



D 



exp(-|x|/LD). 



(11) 
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FIG. 1. Spin polarization drag induced by electric current 
in ID geometry, Eq. (^. The dashed line corresponds to the 
spin polarization profile in the absence of current. The indi- 
cated strength of the electric field, E, is in V/cm. Note that 
the stronger the applied electric field is, the farther away the 
spin polarization can be dragged. The parameters are char- 
acteristic of 10^*^ cm"^ Si-doped GaAs: /i = 3000 cmVVs, 
Ti = 100 ns, T — 1.6 K. Carriers are assumed to be 
non-degenerate. 
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The series of solutions for various strengths of the elec- 
tric field are shown in Figure |^. The parameters are 
characteristic of 10^^ cm"^ Si-doped GaAs ^ — 
3000 cmVVs, Ti = 100 ns, T = 1.6 K. Th e drift cond i- 
tion for this system is realized for E > l^kBT j qpT\ ^ 
2 V/cm. Experimentally, much larger electric fields have 
been applied to this system 

Under spin-drift conditions, the characteristic length of 
the exponential decay of spin polarization is proportional 
to the longitudinal spin relaxation time, Le = VdTi. In 
this regime, there is almost one-to-one correspondence 
between the distance from the interface and the time 
spent by electrons in the region with zero field, x = Vdt, 
which is only weakly altered by diffusion. Therefore, by 
measuring the decay length and knowing the drift veloc- 
ity, Vd = ^lE^ one can directly determine Ti. An ad- 
vantage compared to the conventional pump-probe ap- 
proach H] is that the measurement can be conducted in 
a dc setup that doesn't require ultrafast optics and elec- 
tronics. The spin polarization detection with sub- micron 
spatial resolution can be done by standard optical tech- 
niques based on Kerr/Faraday rotation, or possibly us- 
ing magnetic resonance force microscopy (MRFM) Hby 
modulating the current at the cantilever frequency |10|] . 

In addition to the longitudinal relaxation time, Ti, 
which defines how fast the longitudinal spin polariza- 
tion disappears, an almost identical approach can be ap- 
plied to determine the transverse spin relaxation time, 
T2, which measures the electron spin decoherence time. 
For that, in the setup described above, a weak uniform 
magnetic field, H^, should be applied perpendicular to 
the direction of the strong field that exists at x < 0. In 
the region with the strong field, this additional compo- 
nent will only weakly alter the the direction of the total 
field. However, in the region where originally there was 
no magnetic field, this component will cause the preces- 
sion of the drifting spin polarization. The spatial period 
of the precession 



Lp = 2TTVdluJL-, 



(12) 



is related to the Larmor frequency, oj^ — giisH^/h, 
which is about 62 MHz at Hj_ = 100 G for g = -0.44. 
For the driving electric field E = 10 V/cm and the above 
parameters, Lp 5 fim. Such spatial modulation can be 
resolved both by optical means or with MRFM. To de- 
termine the transverse spin relaxation time T2 one needs 
to measure spatial dependence of the spin polarization 
component perpendicular to H± and then fit it to 



S oc e:sq){—x/vdT2) cos{2Trx / L p) . 



(13) 



Here we neglected the diffusive spread of the spin po- 
larization, which will case an analogue of NMR inho- 
mogeneous broadening. The diffusive spread will matter, 
however, only when it becomes comparable to the preces- 
sion length, Lp. This effectively places an upper limit on 
the distance from the interface for which the precession- 
induced oscillation can be resolved, L < \fiEL'j,/D\ ^ 
250 /im for the above parameters. The upper limit can 



be modified by the choice of the transverse magnetic field 
B± with the only constraint that ujlT2 > 1. Finally, the 
carrier g-factor itself can be determined from the Larmor 
frequency, which can be extracted from the precession 
length. 
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FIG. 2. Spin species separation by inhomogeneous mag- 
netic field. An electric current flows in 2D (film) through a re- 
gion of localized magnetic field, Vo-(r) = — O.Gcr exp[— (r/5)^°], 
with r in micrometers and energy Va in Kelvin. The spin po- 
larization is measured in the units of the total density. The 
strength of the electric field, E, is in V/cm. Notice the re- 
gion of negative ( "minority" ) spin polarization that first forms 
were the current impinges on the region of the magnetic field. 
For large enough drift velocities, the "majority" and the "mi- 
nority" carriers forms a three-tail structure that extends over 
the drift length, Le, downstream. The material parameters 
are the same as in Fig. |l| 



Another application that is only possible in the spin- 
drift regime is the spatial separation of spin species by 
driving electric current through a region with an inho- 
mogeneous magnetic field. To demonstrate this effect, we 
consider a two-dimensional system with a compact region 
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containing magnetic field. Unlike the one-dimensional 
case considered above, the "minority" (opposite to mag- 
netic field) spin carriers can avoid this region by flowing 
around it. Indeed, this is what we observe by solving 
numerically the 2D spin-relaxation drift-diffusion equa- 
tions. Fig. 1^. The calculation is performed on a 100 x 150 
lattice that spans 30 x 45 /im region for parameters rel- 
evant to Si-doped GaAs (specified above). The field 
is localized in 10 /im disk. For large enough electric 
field, there is formation of the negative spin polariza- 
tion where the current impinges on the boundary of the 
magnetic field region [ p4| . The minority spin polariza- 
tion is then dragged downstream by the current until it 
relaxes and diffuses away. Notice that in the zero-current 
state there is no negative spin polarization in the system. 
The "majority" spin polarization is ripped off the local 
field region, and together with the "minority" carriers 
forms a three-tail structure that extends over the spin- 
drift length, Le, downstream. Therefore, driving current 
through this system leads to an effective separation of the 
two spin species. A similar effect can be expected also 
in the ballistic regime realized in the ultra-high mobility 
2DEG structures. 

The effect that we discussed above is analogous to the 
classic Stern-Gerlach effect; however, unlike the original 
effect the spin separation can be achieved even with elec- 
trons. In the original Stern-Gerlach experiment, neu- 
tral spin- 1/2 particles passing through a region with a 
magnetic field gradient get deflected by the force F = 
g^J.B{<J ■ V)H, and produce two spots on the screen that 
represent two possible projections of spin a. For charged 
particles, the situation is complicated due to the Lorentz 
force that acts to scramble the separation. It was first 
pointed out by Bohr that in the case of electrons, the ob- 
servation of Stern-Gerlach effect is impossible because the 
defiection by the Lorentz force exceeds the spin-related 
splitting |l^. The Bohr argument does not apply, how- 
ever, to constrained geometries, e.g. a thin film sample, 
where the Lorentz force is compensated by the Hall volt- 
age. This makes the solid-state Stern-Gerlach effect pos- 
sible. It is important to stress, however, that despite the 
appearance, the mechanism for spin separation here is 
distinctly different from the the classic Stern-Gerlach ef- 
fect. A time- dependent Stern-Gerlach effect was recently 
considered by Fabian and Das Sarma |l^ . 

In the above discussion we assumed a possibility of a 
magnetic field profile that is inhomogeneous but always 
pointing along the same line, e.g. perpendicular to the 
sample plane. This is an important assumption because 
the fringe fields have ability to destroy spin polarization 
through precession. For the 10^^ Si:GaAs, a 100 Gauss 
fringe field at 10 V/cm drive could scramble the polar- 
ization after 5 /im. However, since the fringe fields decay 
rapidly in space, their detrimental effects can be reduced 
by increasing the drive or by reducing the size of the mi- 
cromagnets. Moreover, for thin-film samples, the effects 
of the fringe fields can be essentially eliminated by plac- 
ing the sample between two closely spaced magnet ter- 



minals (split ring) or between two "identical" magnets. 
Alternatively, one could selectively boost the g-factor, 
e.g. by locally doping II- VI semiconductors with Mn. 
This way, when placed even in a relatively weak uniform 
external magnetic field, different parts of the sample will 
feel very different magnetic potentials. By combining 
these approaches, one can engineer a variety of strong 
magnetic potentials which are not affected by the fringe 
fields. 

In summary, we explored the spin-drift regime, which 
is realized for sufficiently large but realistic current den- 
sities in semiconductors. In this regime, the spin density 
can be carried by the current over distances controlled by 
the spin-drift length, which can significantly exceed the 
spin-diffusion length. This opens a possibility of design- 
ing spintronic devices with the spin properties dynami- 
cally controllable by the applied bias. This is in sharp 
contrast with the conventional particle-hole semiconduc- 
tor electronics where electrostatic screening practically 
precludes large scale dynamic manipulation of the rel- 
ative particle-hole density. As specific applications, we 
proposed measurement schemes for the carrier spin re- 
laxation times, Ti and T2, as well as demonstrated the 
possibility of spin species separation by driving current 
through a region with an inhomogeneous magnetic poten- 
tial. The latter is analogous to the classic Stern-Gerlach 
effect. 

Upon completion of this paper, author became aware of 
the related one-dimensional study of the spin-drift regime 
by Yu and Flatte (l|. 

We thank S.A. Crooker, J. Albrecht, D.L. Smith, J. 
Fabian, and S.A. Trugman for useful discussions. This 
work was supported by DARPA SPINs program. 
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